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. Abstract. The purpose of this article is to establish upper and lower estimates for the 

. ' integral kernel of the semigroup exp(— tP) associated to a classical, strongly elliptic pseu- 

1-^^ , dodifferential operator P of positive order on a closed manifold. The Poissonian bounds 

d ' generalize those obtained for perturbations of fractional powers of the Laplacian. In the 

selfadjoint case, extensions to t G C-|- are studied. In particular, our results apply to the 
Dirichlet-to-Neumann semigroup. 
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in ! Introduction 

\q 

^ I Let M be a compact n-dimensional Riemannian manifold and P a classical, strongly 

O ■ elliptic pseudo differential operator (i/jdo) on M of order d > 0. We consider upper and 
lower estimates for the integral kernel }Cv{x, y, t) of the generalized heat semigroup V{t) — 
e~^^ . Semigroups generated by such nonlocal operators have been of recent interest in 
different settings. 

1) For a Riemannian manifold M with boundary M, the Dirichlet-to-Neumann operator 
is a first-order pseudodifferential operator on M with principal symbol |^|. Arendt and 
Mazzeo [AM07], [AM12], initiated the study of the associated semigroup and its relation 
to eigenvalue inequalities, motivating later studies e.g. by Gesztesy and Mitrea [GM09] 
and Safarov [S08]. 

2) The heat kernel generated by fractional powers of the Laplacian A''/^ and their 
perturbations provides another example. Sharp estimates for e~*^ , < (i < 2, can be 
obtained from those for e~*^ by subordination formulas. For perturbations on bounded 
domains in M"^, recent work on estimates includes Chen, Kim and Song [CKS12] and other 
works by these authors, and Bogdan et al. [BGRIO]. 

In this article we generalize the Poissonian estimates obtained in the second case to 
parameter-elliptic operators P on closed manifolds, by pseudodifferential methods. In 
particular, we allow nonself adjoint operators. A main result is: 
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Theorem. The kernel of the semigroup satisfies 

(*) \K:vix,y,t)\<Ce-'''H{d{x,y) + t^/'^)-''-'^, for x,y e M,t > 0, 

for any ci smaller than the infimum j{P) of the real part of the spectrum of P. 

If P is self adjoint > 0, the estimates extend to complex t = e*^|t| for \9\ < ^, with 
uniform estimates 

(**) \K:v{x,y,t)\ < C(cos^)-^e-^(^)^«* {d{x,y) + ^^^^^'^^ + 1^1'^'^"'' + 

where N = max{|, ^ + 4d + 7}. 

Here d{x, y) denotes the distance between x and y. If P is a system, it suffices that 
P — A is parameter-elliptic on the rays in a sector containing {Re A < 0}. Extending (*), 
also derivatives of the kernel, and, if further spectral information is available, a refined 
description of the long-time behavior, are obtained in the paper. 

For the Dirichlet-to- Neumann operator, as well as for the perturbations of fractional 
powers of the Laplacian of orders < d < 2, we get not only upper estimates but also 
similar lower estimates at small distances. 

The estimate (*) exhibits a large class of operators which satisfy upper estimates closely 
related to those studied abstractly in Duong and Robinson [DR96] . As a simple application 
of (**) and Holder's inequality, one can for example obtain ultracontractive estimates 

uniformly for t G C+. In the case of operators with Gaussian heat kernel estimates, a rich 
spectral theory has been developed (see e.g. Arendt [A04], Ouhabaz [O05]). 

With the help of comparison principles, our result implies Poissonian estimates e.g. for 
boundary problems in an open subset O of M: If P is the variational operator associated to 
a Dirichlet form a with domain V C L2(M), we consider the abstract Dirichlet realization 
Pq, associated to the closure of a|i?nCo(n)- ^^e case where a is Markovian, one obtains 
< JC^-tPn < JCg-tp on Q. See Grigor'yan and Hu [GH08] for more refined comparison 
principles. 

Outline. Section 1 collects some known facts. In Section 2 we treat semigroups generated 
by nonselfadjoint P for t > 0. Section 3 extends the estimates to complex t for selfadjoint P. 
Section 4 includes lower estimates for perturbations of fractional powers of the Laplacian 
and for the Dirichlet-to-Neumann operator. 

1. Preliminaries 

Notation: {$,) = ^/$,^ + 1. The indication < means "< a constant times" , > means 
"> a constant times" , and = means that both hold. 

Let P be a classical ipdo of order d e R+, acting in a Hermitian AT-dimensional C°° 
vector bundle E over a closed, compact Riemannian n-dimensional manifold M. We 
assume that P — A is parameter-elliptic on all rays with argument in ]^ — (po, ^ + (po[ for 
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some (fio e]0,^[. (In the notation of Seeley [S67] , these rays are "rays of minimal growth" .) 
From P one can define the generahzed heat operator V{t) = e~*^, t > 0, a holomorphic 
semigroup generated by P, as explained in detail e.g. in [G96], Sect. 4.2. The kernel 
ICv{x,y,t) (C°° for t > 0) was analyzed there in its dependence on t, but mainly with a 
view to sup-norm estimates over all x, y, allowing an analysis of the diagonal behavior, that 
of )Cv{x, x,t). We shall expand the analysis here to give more information on /Cy (x, y,t). 

For convenience of the reader we recall the definitions of symbol spaces that are used. 
For deR, the symbol space Slo{W x R") consists of the C°° -functions a{x, C) {x, C e M'^) 
such that for aU Q!,;5 e Ng', 

(1.1) i^fa(^, 01 <(0 

it is a Frechet space provided with the seminorms sup^^^ \{^)~^'^\°'\D^D^a\. The symbols 
define operators A = Op{a{x, ^)) of order d by 

Op{a{x,0)u = :F-^a{x,OJ'u= [ e-^a(a;, 0^(0 ^C, 

where J^u = u denotes the Fourier transform and — (27r)~"(i^. The operator maps 
from 5(11^"^) to 5(M"^), extending to suitable spaces of distributions and Sobolev spaces, 
and obeying various composition rules. 

The space of classical symbols of order d, S'^{W x M"), is the subset of Sf^Q{W x W^) 
where a{x,^) moreover has an asymptotic expansion a ~ X^^gj^g aa-i in terms ad-i{x,^) 
homogeneous in ^ of degree d — / for |^| > 1, such that a'j^ = a — ^i^m ^d-i £ Sf^^ for 
all M e Nq. The principal symbol is often denoted aP . 

It should be noted that we here use the globally estimated symbols of Hormander [H83] , 
Section 18.1, which have the advantage that remainders are kept inside the calculus. 

Operators on manifolds are defined by use of local coordinates and rules for change of 
variables, composition with cut-off functions etc.; we refer to the quoted works for details. 

The book [G96] moreover includes parameter-dependent symbols a(x, ^, A) for A in a 
sector of C, with special symbol estimates involving the parameter (also operators on 
manifolds with boundary are treated there). 

Consider a localized situation where the symbol p{x,^) of P is defined in a bounded 
open subset of M"^ — we can assume it is extended to M*^, with symbol estimates valid 
uniformly in x. The abovementioncd hypothesis of parameter-ellipticity means that the 
spectrum of the principal symbol p^\x,^) (an x A^-matrix) is contained in the sector 
{A I I arg A| < 6*0}, 6*0 = f — (/?o, when |^| > 1. This holds in particular when P is strongly 
elliptic, for then 

(1.2) Re(j9°(a;, C)v, v) > c|^|'^|^;p, for |^| > 1, w e C^, with c> 0, 
and hence since 

(1.3) |Im(A,^^)| < |(A,^^)| < C\i\'^\v\^ < c-^CRe{p^v,v), for |^| > 1, ^; G C^, 

P satisfies the condition of parameter-ellipticity with (/jq = §^^07 where ^0 = arctan(c~^C) 
e]0, ^[. When P is scalar, the two ellipticity properties are equivalent, but for systems. 
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strong ellipticity is more restrictive than the mentioned parameter-eUipticity (also called 
parabolicity of dt + P). 

The spectrum (t{P) of P lies in a right half-plane and has a finite lower bound 7(P) = 
inf{ReA | A G a{P)}. We can modify for small ^ such that cr(p°(a;,^)) has a positive 
lower bound throughout and lies in {A = re*^ | r > 0, 16*1 < 6*0}. 

The information in the following is taken from [G96], Section 3.3. 

The resolvent Q\ = {P — \)~^ exists and is holomorphic in A on a neighborhood of a 

set 

(1.4) Wr^,e = {\&C \ |A| >ro, argAG [^0 + £, tt - 6*0 - e], ReA<7(P)-£}. 

(with £ > 0). There exists a parametrix on a neighborhood of a possibly larger set 
(with 5 > 0, £ > 0) 

V5,e = {A G C I |A| > 5 or arg A G [6*0 + £,7r - 6*0 - e]}; 

such that this parametrix coincides with (P — A)~^ on the intersection. Its symbol g(x, ^, A) 
in local coordinates is holomorphic in A there and has the form 

(1.5) q{x,^,X) r^^q-d-l{.x,^,X), where = (/(x, - A)"^ 
Here when P is scalar, 

21 

(1-6) = 6i,i(a;,C)?id, = ^6«,fe(a;,09^d\ ••• ; 

with symbols hi^k independent of A and homogeneous of degree dk — I m. ^ for |^| > 1. 
When P is a system, each q-d-i is for Z > 1 a finite sum of terms with the structure 

(1-7) r(x, A) = h^q''_\h2q''_\ ■ ■ ■ 6m?::^6m+i, 

where the are homogeneous tpdo symbols of order Sk independent of A, the are positive 
integers with sum > 2, and • --I-sm+i— <^(i^i + - • •+i^M) = —d—l. (Further information 
and references in Remark 3.3.7.) Moreover, the remainder q'^ — q — ^i<M Q-d-i satisfies 
for A with |7r — arg A| < ^ + (fi, any < (po, 

(1.8) \DlD^q'M{x,C,X)\ < + + |A|Vd)-2d^ ^hen M + |a| > d. 

(Cf. Theorems 3.3.2, and 3.3.5, applied to the rays with arguments in ]^ — ipo, ^ + (Pq[.) 

2. Semigroups generated by parameter- elliptic 
pseudodifferential operators 

As explained in [G96], Section 4.2, the semigroup V{t) = e~^^ can be defined from P 
by the Cauchy integral formula 

(2.1) v{t) = ^ j^e-'\P-\)-U\ 
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where C is a suitable curve going in the positive direction around the spectrum of P; it can 
be taken as the boundary of Wro,e for a smaU e. In the local coordinate patch the symbol 
is (for any M e No) 

v{x,^,t) = v-d H \-V-d-M+i +v'm ^ ^v-d-i{x,^,t), where 

(2.2) 



v-d-l = ^ / e ^^q_d-i{x,^,X)dX, v'm = ^ e '^qudX. 
Jc Jc 



A prominent example is e~*^^ where A denotes the (nonnegative) Laplace-Beltrami 
operator on M. This is a Poisson operator from M to M x as defined in the Boutet de 
Monvel calculus ([B71], cf. also [G96]), when t is identified with x^+i. When M is replaced 
by R", its kernel is the well-known Poisson kernel 

(2.3) ]C{x,y,t) = Cn- 



{x-y,t)\^+^ 

for the operator solving the Dirichlet problem for A on 

Also more general operator families V{t) = e~*^ with P of order 1 are sometimes 
spoken of as Poisson operators (e.g. by Taylor [T81]), and indeed we can show that for 
P of any order d e M_|_, V{t) identifies with a Poisson operator in the Boutet de Monvel 
calculus. This will be accounted for in detail elsewhere. In order to match the conventions 
for Poisson symbol-kernels, the indexation in (2.2) is chosen slightly differently from that 
in [G96], Section 4.2, where V-d-i would be denoted V-i. We define V-d-i{t) and Vl^{t) in 
local coordinates to be the i/jdo^s with symbol v-d-i{x, ^^t) resp. t'^(x,^,t). The kernel 
ICv{x,y,t) is in local coordinates expanded according to the symbol expansion: 

(2.4) Kv{x,y,t)= )Cv_a-i{x,y,t) + K:v;^ix,y,t). 

0<1<M 

The following result follows from [G96]. 
Theorem 2.1. 1° In local coordinates, the kernel terms satisfy for some c' > 0: 

f{l-n)/d 

(2.5) \lCv_,_Ax,y,t)\<e-'''{ t{\\ogt\ + l)ifd-l = -n, 

t if d — l < —n. 

For a given cq > we can modify to satisfy infa;^^ 7(^*^(0;, ^)) > cq; then c' can be any 
number in ]0, co[ . 

2° Moreover, with the modification in 1° used with cq — j{P) ifj{P) > 0, the remainder 
satisfies 

t(M-n)/d ifd-M>-n, 

(2.6) \]Cvi,{x,y,t)\ ti\logt\ + l) ifd-M = -n, 

tifd-M<-n, 
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for any c\ < 7(-P). In particular, 

(2.7) \)Cv{x,y,t)\<e--^H-^/''. 

Proof. The theorem was shown with shghtly less precision on the constants c', ci in [G96], 
Theorems 4.2.2 and 4.2.5. It was there aimed towards apphcations where d is integer. 
The estimates of resolvent symbols in Section 3.3 are still valid when d e , but the 
replacement of P by P + a (a G M) in the beginning of Section 4.2 on heat operators only 
gives a classical i^xlo when d is integer, so we need another device to take the value of j{P) 
into account for general d e M+. We shall now explain the needed modifications, with 
reference to [G96]. 

For 1°, the proof in Theorem 4.2.2 shows the validity of (2.5) with a small positive 
c' < infa;^^ 7(p*^(x, ^)). For a given cq > 0, the proof goes through to allow any c' < cq, 
when p^{x, ^) is modified for |^| < R (for a possibly large R) to satisfy inf '^{p^{x,^)) > cq. 

For 2°, the remainder symbol q'j^ is holomorphic on Wro,s', here if 7(-P) > we define 
the terms q-d-i as under 1°, with cq = 7(-P)- For large M, q'j^ is < (A)~^. The proof of 
Th. 4.2.2 gives an estimate of /Cy^ by e~'^i*t (1 + | logt|), and the proof of Theorem 4.2.5 
shows how to remove the logarithm. The estimates of ICy' for lower values of M follow 
by addition of the estimates of finitely many JCv-d-i'^^^^^- ^ 



We shall improve this to give information on the dependence on |x — y| also. This will 
rely on the following result on kernels of SI Q-ipdo^s, found e.g. in Taylor [T81], Lemma 
XII 3.1, or [T96], Proposition VII 2.2. 

Proposition 2.2. Let a E C°°(M" x W) be such that for some r e R, N e with 
N >n + r, and all < |a| < N, 

(2.8) sup,,^(0-^+l"l|i?fa(x,OI<oo. 

Then the inverse Fourier transform K.A{x,y) = J-^^^a{x,^)\z^x-y is 0(\x — y\~''^) for 
\x — y\ ^ oo, and satisfies for \x — y\ > 0; 

\x — y\~'^~^ if f > —n, 

(2.9) \JCA{x,y)\<{ |log|a;-y|| + 1 i/r = -n, 

1 if r < —n. 

In particular, if a E SIq{M.'^ x M"') defining the ipdo A, the estimates hold for its kernel 
)Ca{x, y) for all N > n + r, each estimate depending only on the listed symbol seminorms. 

The dependence on N follows from an inspection of the proof. 

In the scalar case the kernel study can be based on nice explicit formulas, that we 
think are worth explaining. Consider the contribution from one of the terms in (1.6). As 
integration curve we can here use Cs consisting of the two rays re*^ and re~*^, 9 = 9q + e. 
For i > 0, a replacement of tX by g gives: 



(2.10) 



(^<:f\- ^ ! ^-tA bi,k{x,i) , [ t%,k , 

C — Q 

= 4-t\ k / , / dg = U'^bi ke-*P° ; 
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here we have replaced the integration curve by a closed curve Cg^R connecting the two 
rays by a circular piece in the right half- plane with radius R > 2t\p^{x,^)\, and applied 
the Cauchy integral formula for derivatives of holomorphic functions. This shows: 

21 

(2.11) v-d = e-'P\ v-d-lix,^,t) = J2 M*'kfe(^>e)e-'^°^"'«) for I > 1. 

fc=i 

Then the kernels of the V-d-i (t) can be estimated by the following observations. 

Proposition 2.3. Letp^{x,^) be the principal symbol of a classical scalar strongly elliptic 
i(jdo P on of order d e IR+, chosen such that Rep^{x,^) > cq > 0. 

1° Let c' G [0,co[. For any j e No, (t(p°(a;,0 - c'))^e-*(P°(^'^)-'^') is in Slo{W x W) 
uniformly in t >0. 

2° Let 

(2.12) -^''^'*^ = ^'Ja/' \»M%-^ ''' 
where k > 1 and b E Sf^^'^W x W). Then 

(2.13) w{x,^,t) = ^t'=6(a;,0e"*P"^^'^) = e-'''Hw'{x,^,t), 



—1 



where w'{x,^,t) E Sf^\W x W), uniformly for t > 0. 
Moreover, w{x, z, t) = J-^^^w satisfies for any c' E ]0, co[ 



{t l^l^-^-'^ ifd-l>-n, 
t (I log 1^1 1 + 1) ifd-l = -n, 
t if d — I < —n. 

It follows that for I > 1, lCv_^_i{x,y,t) = T'^^^V-d-l{xi^it)\z=x-y satisfies the estimates 

t\x — y\'--^-'^ if d — I > 

(2.15) \]Cv_,_, {x, y,t)\< e-'''' \ t {\log\x - y\\ + 1) if d - I = -n, 

t if d — I < —n. 

Proof. 1°. For each fixed t > 0, e~*^ ^^'^'^ is rapidly decreasing in ^, hence is in Si^. But 
for our purposes we need estimates that hold uniformly in t for t — )■ 0. Let 

Mj,k,l = sup s^d^is'e-'). 

s>0 

Then for t > 0, ^ e M", 

\{tp'{x,Oye-'P°(-'^^\<M,,o,o, 
(2-16) |%((t/)^e-*P")| = \ds{s^e--^)U^,,otd^y\<M,,kA\{pT'd^.p'\ < iO'^ ■■■ 
\d^{{tpye-'P°)\<{0-^"\,... 
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showing the assertion for c' = 0. (2.16) holds also if p° is replaced by p° — c' throughout, 
when c' e ]0, co[ . 

2°. The first identity in (2.13) was shown in (2.10). We can also write 

w{x,^,t) = it6(/ - c')'-^t(/ -c'))'=-'e-='*e-*(^'°-'=') = e-'''Hw'{x,^,t). 

Here 6(p° - c')^''' is in Sf'^^ independent of t, and by 1°, (t(/ - c'))''-^e-*(P°-'^') is 
uniformly in S^ q, so it follows that w' is uniformly in Sf^K We can now apply Proposition 
2.2 to draw the conclusion (2.14). 

Since V-d-i{x, ^, t) is a sum of such terms when / > 1, the estimates (2.15) follow. □ 

For systems P we can use systematic estimates from [G96] . We find for general P: 

Theorem 2.4. 1° In local coordinates, /Cv_^ satisfies for some d > 0; 

(2.17) \K:v_,{x,y,t)\<e-'''H\x-y\-''-^. 

For I > 1, the kernels fCv-d-i satisfy (2.15). If ^(P) > 0, we modify to satisfy 
infj;^^ 7(p°(a;, ^)) > 7(-P), then c' can he any number in ]0,7(P)[. 
2° Moreover, with chosen as in 1° , 

{t\x - y\^-'^-'' ifd-M>-n, 
t{\\og\x-y\\ + l) ifd-M = -n, 
tifd-M<-n, 

for any ci < 7(-P). In particular, 

(2.19) \K:v{x,y,t)\ < e-''^H\x-y\-'^-''. 



Proof. 1°. When P is scalar, the estimates in (2.15) for / > 1 are shown in Proposition 2.3, 
when we take cq = 7(-P) if l{P) > 0. For general systems P, the symbols q-d-i are sums 
of symbols as in (1.7), and we apply [G96], Lemma 4.2.3. Here (4.2.35) with k = —d — I 
shows that 

for all a, (3. Actually, the estimate (4.2.35) has e~^*^^^ with a positive c as the last factor, 
but an inspection of the proof (the location of integral contours) shows that e~'^*^^^ can be 
replaced by e"*^'*, if c' < inf ■^{p^{x, $,)). This shows that e^'*t~^V-d-i is in ^^q^ uniformly 
in t, so the estimates of the }Cv_a-i follow by use of Proposition 2.2. 

For Z = 0, we can argue as follows in the scalar case: For each j = 1, . . . , n, 

d^.v.d = %e-*^° = -t{d^^py-'p\ 

where d^^p^ e 'S'i,o^- Now as in Proposition 2.3, e~'^'^d^.p^e~*^P°~'^'^ is in Sf^^ uniformly 
in t, and hence V-d = J-^^^V-d satisfies, since d — 1 > — n, 

(2.20) \zjV-d\ < e-"'*t|^r^+^-". 
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Taking the square root of the sum of squares for j = 1, . . . , n, we find after division by 1^1 
that 

(2.21) <e~'''H\z\-'^-''. 
In the systems case we note that 

(2.22) d^.q-d = -q-d{d^,p'')q-d, 

since d^.[{p^ — X){p^ ~ ^)~^] =0- Lemma 4.2.3 apphes to this in the same way as above, 
showing that 

\D^,D^d^^v_dix,U)\ < (O'^-^-l^lte-'^'*, 

so e'^^td^.V-d is is uniformly in S^^^. We conclude (2.20), from which (2.21) follows, 
implying (2.17). 

2°. Here the estimate in (2.18) has already been shown for large M in Theorem 2.1. For 
lower values of M, we can add the estimates of the entering homogeneous terms lCv_^_i 
with / > M; the top term gives the weakest estimate. (It is used that x and y need only 
run in a bounded set, for the contribution from the localized piece.) □ 

Theorems 2.1 and 2.4 together lead to Poisson-like kernel estimates: 

Theorem 2.5. 1° One has in local coordinates: 

t {\x - y| + t^/dy-d-n lffl_ly_^^ 

(2.23) \lCv_,_, {x, y,t)\< e''"'' { t{\ \og{\x -y\+ tVrf)| ^i)ifd-l = -n, 

t if d — l < —n, 

for some c' > 0. Ifj{P) > 0, we modify to satisfy mix, ^jip^{x,^)) > 7(-P); then c' can 
be any number in ]0,7(P)[ . 

2° Moreover, with p'^ chosen as in 1° , 

{t{\x-y\+t^/'^)^-'^-'' ifd-M>-n, 
t (I log(|x -y\+ ti/d)| + l)ifd-M=-n, 
t ifd-M<-n, 

for any c\ < 7(-P). In particular, 

llCvix, y, t) I < e-'^H {\x -y\+ f/d)-d-n^ 
llCv^ix, y, t)\ < e-'-H {\x-y\+ t^/df-d-n^ 



(2.25) 



3° For the operators defined on M, one has (with d{x,y) denoting the distance between 
X and y ) 

(2.26) \lCv{x,y,t)\ < e-^^H{d{x,y) + ^ , 

for any ci < 7(-P). 
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Proof. l°-2°. In the region where \x — y\ >t^/'^, 

\x-y\<\x-y\ + t^^'^ < 2\x - y\, 

in other words, |x — y| = \x — y\ + t^/^. Then the estimates in Theorem 2.4 imply the 
validity of the above estimates on this region. 

In the region where |x — |/| < t^^'^, we have instead that t^^'^ = |a; — |/| + t^/*^. Then the 
estimates in Theorem 2.1 imply the above estimates on that region; for example 

= ^ ^ ^ (l^. _ ^1 + ^1/dyd-n 

there. For the two regions together, this shows (2.23)-(2.25). 
3°. This follows from the estimates in local coordinates. □ 

When the eigenvalues of P with real part equal to 7(-P) (necessarily finitely many) 
are semisimple (i.e., the algebraic multiplicity equals the geometric multiplicity), we can 
sharpen the information on the behavior for t — )■ oo: 

Corollary 2.6. Assume that all eigenvalues of P with real part ^{P) are semisimple (it 
holds in particular when P is self adjoint) . Then 

(2.27) |/Ce-..(-, y. t)\ < --'^''^ \^^,J^^ydY V) + t'^r- + l) • 

Proof. The spectral projections 11^ = j^XP — \)~^d\ onto the eigenspaces Xj for 
the eigenvalues {Ai,...,Afc} with real part 7(-P) (where Cj is a small circle around the 
eigenvalue), are pseudodifferential operators of order — oo, and their kernels ]Cujix,y) are 

bounded. If e > 0, the operator P' = P + e Y^'^^i Hj satisfies 7(-P') > 7(-P)- By Theorem 
2.5 applied to P', 

|/C,_.P. {x, y, t) I < e-^(^)*t((i(x, y) + tVd)-d-n. 



On the other hand, V{t) = e"*^' + (1 - e"^*) ^Li e'^^m 



k 



so 



From 



/Ce-tp (x, y, t) = /C,_tP' (x, y, t) + (1 - e-'') ^ e"*^^ /Cn, {x, y) . 



1 - e-^* < min{l, et} < — — ^ -j-.-: < ^ 



(diam(M) + tVd)d - {d{x,y) + t^/d)d' 

we conclude that (1 - e-^*)|/Cn^.(a;,?/)| < ^a{x,y)+t^/dy ^ ^^'^ (2-27) follows since |e-*^^ | 
e-*7(-P) for each j. □ 
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Remark 2.7. The proof of Corollary 2.6 allows to sharpen the estimates in Theorem 2.5 
and Theorem 2.9 below also in the general case where the eigenvalues with real part 7(-P) 
are not all semisimple. Denote by r the dimension of the largest irreducible P-invariant 
subspace of any eigenspace Xj associated to an eigenvalue with real part 7(-P). Then in 
Theorems 2.5 and 2.9 we may replace the upper bound e-'''H{d{x,y) + t^/'^)-'^-''-^ by 

(2.28) e-(^)*(l + t'-') ^d^,J^^^/,y m-^v) + t'^')—' + l)- 

It is not hard to extend the estimates to complex t in a sector around 'R+ . Namely, since 
p° has its spectrum in the sector {| argA| < ^o}, e*'^P satisfies the parameter-ellipticity 
condition when \(p\ < (pQ = ^ — Oq. For each (f it generates a semigroup e~*^ and these 
operator families coincide with the holomorphic extension of V{t) to the rays {re^'^} in the 
sector = {t G C | | argt| < Lpo}. On each ray we have the estimates in Theorem 2.5, 
they hold uniformly in closed subsectors of V^g . We have hereby obtained: 

Theorem 2.8. With ^po and 9o defined as in the beginning of Section 1, the semigroup 
generated by P extends holomorphically to the sector {| argt| < ipo), and the estimates in 
Theorem 2.5 hold in terms of \t\ on any closed sector {|argt| < p} with < <^ < (fQ, 
taking ci < min|^/|<(p 7(6*''' P). 

More information in the case where P is selfadjoint will be given in Section 3 below. 
Also the derivatives of the kernels can be estimated by use of the symbol estimates in 
[G96]. 

Theorem 2.9. 1° One has in local coordinates: 

(2.29) \D^,D^yDiJCv_,_,{x,y,t)\< 

{t {\x -y\ + t^/dy-a+j)d-h\-n (^j ^ + _ ; > _n, 
t (I log(|a; -y\+ t^^) \ + l)if {j + l)d + I7I - ^ = -n, 
tifij + l)d+\^\-l<-n, 

for some c' > 0. If'y{P) > 0, we modify p'^ to satisfy infa;^^ 7(^^(0;, ^)) > 7(-P); then d can 
be any number in ]0,7(P)[. 

2° Moreover, with p^ chosen as in 1°, 

(2.30) \DPD^yDllCv;^ix,y,t)\< 

( t{\x-y\+ tl/d)M-0-+l)d-|7|-n ^ ;l)c/ + |7| - M > -U, 

e-^i* I t{\ log{\x -y\+ t^/^)! + 1) if {j + l)d + I7I - M = -n, 
[ t if {j + l)d + I7I - M < -n, 

for any ci < 7(-P). 

3° The estimates of derivatives o//Cy hold for the operator defined on M with \x — y\ 
replaced by d{x,y). 
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Proof. As in Theorem 2.7, the estimates are pieced together from estimates generahzing 
those in Theorem 2.1 resp. Theorem 2.4 to include derivatives. We use that 

To generahze Theorem 2.1 to allow x- and y-derivatives we just have to apply the 
arguments of [G96] Theorems 4.2.2 and 4.2.5 to the modified symbols ^'^D^v-d-i, to 
get the estimates (2.29) with \x — y\ replaced by 0. Derivatives with respect to t alone 
are explained in Theorem 4.2.5; finally this is combined with x- and y-derivatives in a 
straightforward way. Similar considerations work for remainders; here we can in fact refer 
directly to (4.2.60) for large M, and the statements for lower M follow by addition of the 
appropriate set of estimates of }Cv_a_i-terms. This gives the expected generalization of 
Theorem 2.1, namely (2.29)-(2.30) with \x — y\ replaced by 0. 

For the generalization of Theorem 2.4 we note that estimates 

\CD^,D^Div^d-lix,^,t)\ < ^^)U+l)d+h\-\a\-l^^-c't 

for |q!| +Z > 0, all /3, j, follow from [G96] Lemma 4.2.3 (sec the remarks around (4.2.40) for 
how to include t-derivatives, as done also in Theorem 4.2.5). Thus e'^ ^t~^^'^' D^Dlv-d-i is 
in uniformly in t, and it follows by Proposition 2.2 that 

r t\z\^(j + ^)d-\l\+l-n^ if + _/ > _n, 

\D2D^^Dlv-d-l{x,z,t)\ < e-'U t (| log(|^| + tV^)| + i) if + i)d+ h\ - I = -n, 

[ t if {j + l)d+ |7| - / < -n. 

This implies estimates as in (2.29) with \x — y \ + 1^^^ replaced by \x — y\. The conclusion 
is immediate for Z > 1, and for Z = 0, we use the estimates of D^.v as in the proof of 
Theorem 2.4. Again for remainder estimates, we can appeal to (4.2.60) for large M. 
The proof is now completed as in Theorem 2.7. □ 

3. Estimates in the complex plane for selfad joint operators 

In this section we shall derive some uniform estimates for the extension of the semigroup 
into the region C_|_ = {t e C | Ret > 0}, when P is selfadjoint. To do so, we need to 
account for how the estimates of symbols and remainders like (1.6-8) depend on arg A. We 
assume for simplicity that P is > 0. 

As already observed in Theorem 2.7, V{t) exists for t G C+. There are uniform estimates 
on closed subsectors, but to describe the behavior for rays near the imaginary axis we need 
estimates of q{x,^, A) for A near R+. 

As in [G96], we denote |A|^/'^ = and write ((C,A*)) = (1 + ICP + A*^)^^^ for short as 
(^,/.);itis =(l + |^| + |A|Vd). 

Proposition 3.1. Let P be selfadjoint > and let X e C with arg X = (p e]0,^[. Then 

\q-d{x,C, A)| = \{p'{x,0 - X)-'\ < (sin^)-i(^,/.)-^ 
|Dfi^fg_d(x,e,A)|<(sin^)-i-H-l/3|(^)d-'-M(^,^)-2d, when \a\ + \P\ > 0. 
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For all I, a, j3 with I > 0, 

(3.2) |Dfi?fg_d-K^,^,A)| < (sinv:)-2«-l«l-l/?l(^)'^-'-l«l(^,^)-2<i. 

Proof. We have for A = e*'^|A| with (p e]0, |-[, v e C"'^, since p^{x,^) is symmetric with 
lower bound > c(^)'^: 

\{p^v, v) — X\v\'^\ > I Im((p°f , v) — \X\e^'^\v\'^)\ = |A| sin(/7|u| 



\v\ 


|2| > 1 


\v\ 


2| _ 



= sm(p{p^v,v) > si-n(pc{^)'^\vf, 

from which follows 

- A)^||^| > |((/ - A)^,^)| > sin^(|A| + iO'M'. 

This implies that |(p°-A)-i| < {sm(fi)-^{{C)'^ + \X\)-^ = {sm(fi)-^{^, /i)-'^, showing (3.1). 
The other estimates follow as in [G96] from the structure of the terms in the parametrix, 

using (3.1): q-d-i is for Z > 1 a finite sum of terms, where ui -\ h um > 2 takes values 

up to 21, 

r{x,t A) = hq''J^b2qZ ' ■ ■ buQ^'J^bM+u 

of. (1.7). Each Q-d contributes with a factor (sint/?)"^, and there are up to 21 such factors; 
this shows (3.2) for ct = /3 = 0. Each differentiation may hit a factor q-a giving an extra 
(sinf^)"-*^ in view of (2.22); this leads to the estimates (3.2) by the Leibniz formula. □ 

Estimates of remainders Q'j^ and their symbols are more difficult to work out, since 
they depend on the interplay between the exact resolvent Qx and the homogeneous symbol 
terms, and they will be more costly in powers of {sia(p)~^ the larger M is taken. We shall 
here go directly to remainder kernel estimates. 

More precisely, we consider the kernel of the operator Q'j^ — Qx — ^^^^ Q-d-i, where 
each Q-d-i is an operator on the manifold M constructed from the symbols q-d-i in local 
coordinates. 

For A as in Lemma 3.1, 

||(P- A)m||||m|| > |((P- A)^,^)! > ImAll^lP = sinful A| \\uf; 

hence the resolvent Qx — {P — has operator norm < {s'mif \\\)~^ in L2{M). The 
operator Q'j^ = Qx — Y1iI<m Q-d-l is a ipdo of order —d — M (for each A). Moreover, 

Q'm = Q'm{P - ^)Q\ = RmQx, where 
^^■^^ Rm = Q'm{P-X)^1- J2^^^Q-d-i{P - A) 

is a Ipdo of order —M. Now 

\\RmQx\\c(L2,H^) < (sinv? |A|)~"^||i?M||z:(L2,H«)- 

For s > n, it is known that a ijjdo continuous from L2(M) to H^{M) has a continuous 
kernel estimated by the operator norm; then the kernel of Q'j^ = RmQx is continuous and 
is estimated by 

(3.4) \ICQ>Jx,y,X)\ < \\Q'M\\ciL„Hs) < {smip\X\)-'\\RM\\ciL„H^)- 

In preparation for the study of Rm and its dependence on A we prove a lemma on a 
typical composition formula, treated by basic methods in the ipdo theory. 
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Lemma 3.2. Let b{x, e 5fo(M" x R^), and let a{x, A) G S^lo{W x M") with respect 
to (x,^), with A as in Proposition 3.1, such that for some d' > 0, N & M. one has for all 

(3.5) \D^D^a{x,^,X)\ < (sin^)-^-l«l-l^l(0'^'+'^^"l"l(^,/x)-'^' . 

1° There exists c{x,^,X) G ^fo+'^2(M" x W) such that Op(a) Op(6) = Op(c), and for 
every M e No, 

(3.6) c{x,^A) = J2^ Dfa{x,^,X)d^b{x,0 + rM{a,b) , 
where 

(3.7) \D^DfrMia,b)\ < (sin^)-^-^-l"l-l^l (O'^'+'^^+^^-'^-l'^l (e, Z^)"'^'. 
2° /f (3.5) for a = (3 = is replaced by 

(3.8) |a(x,e,A)| < (sinvp)-^(0'+'Me,/^)-' • 

/or some < d < d', then (3.6) /loZds with (3.7) va/zd /or M > 1 and the estimates of ro 
replaced by 

(3.9) |L>fDfro(a,6)| < (siiK/?)-^-!"!-!/^! (O^+'^^+'^^-l"! (e, /z)"^. 
3° For 7 e N^, D'^ Op(a) = Op(a^), where 

(3.10) |i?fi?^"a^(x,e, A)| < V^^, (sin^)-^-'=-l"l-l/5|(e)'^'+'^^+l^|-'=-|«l(^,/.)-'^' . 

'A;<|7| 

Proof. 1°. Let x(x,^) denote a C°°-function that is 1 for \x\'^ + < 1 and vanishes for 
I^P + I^P > 2, then we can replace the given symbols by their products with x(£x,£^), 
which makes all integrals calculated below convergent. It is known in the theory (by the 
technique of oscillatory integrals, cf. [H83], Sect. 7.8), that the resulting symbols converge 
to the given symbols for £ — )■ in all the seminorms that are involved. The modified 
symbols will again be denoted a, b. We can also assume that b has compact support in 
X (in a set containing the x for which we need the formula). Then b{r],^) = J-'x^nb{x, ^) 
satisfies 

(3.11) \Dfb{v,0\<{v)-'''{0''-^''^, 

for all a, N' . It follows from the '0do defining formula that Op(a) Op(6) = Op(c), where 
c{x, A) = [ a{x, r], X)b{y, ^)e^(^-2/)- V(^-^>« dzd^dydrj 

(3.12) 

= / a{x,^ + ri,X)b{ri,Oe''"''dr]. 
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If M > 0, we insert the Taylor expansion of a in ^ up to order M, 

+ E|«|=mS^" l'^^ - h)^-'d^a{x, e + hr), A) dh, 
obtaining that c = c<m + tm, where 

c<M = {27Tr- [ V ^d^a{x,^,X)v''b{v,Oe'^''dv 

The sum over \a\ < M equals the sum in (3.6). For the last integral we use that 

\d^a{x,^ + h7],X)\ < {smip)-^-^'{^ + hr]f+^'-^{C + hri,fj,)-'^' 

< (sin ^)-^-^ (0'^'+^^ (e, 1^)-'' iv) , 

by the Peetre inequality. Taking this together with the estimates (3.11) of b (with a large 
N'), we can conclude that 

\rM\ < (sin<^)-^-^(0'^'+'^^+^^-^(e,/^)-''. 
For M = 0, we apply such considerations directly to ro = c{x, ^, A) in (3.12): 

|ro| < (sin^)-^|(e + r/)^'+'^^(e + r/,/x)-^'(r7)-^'(e)^^^rr7 < (sin<^)-^(e)'^'+'^^+'^Me,/^)-''. 
Derivatives of tm in x and ^ are treated in a similar way. 

In the case 2° the proof goes through in a similar way, except that d' is replaced by d 
in expressions containing undifferentiated factors a. 

In 3°, the A-independent factor is to the left, and (3.12) holds with integrand 
(C + vV ^z-^rfi{z,^,X)e^^''^ . The Taylor expansion of (^ + rij)"' is a finite binomial ex- 
pansion X1k<7 (k)'?^ '^^'^ '^'^'^ leads to a finite composition formula where the estimates 
(3.10) of the terms can be read off directly. □ 

The composed symbol c = ro(a, h) is also denoted a oh (used in [G96]) or a#6. 

The lemma will be used in the following investigation of the symbol of Rm- We de- 
note P — X — P, with the parameter-dependent symbol p{x,$,,X) — p{x,$,) — A in local 
coordinates; here for any M e Nq, 

P = ^ Pd-k+PMi = 5Z Pd-k+PMi with 

(3.13) k<M k<M 

Pd=p- X, pd-k = Pd-k for /c > 0, p'm = p'm Af > 0. 

Pd is also denoted p^ . The Pd-k are homogeneous in |^| of degree d — k iov \^\ > 1, and 
p'm e Sf-Q^{W xW). 
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Proposition 3.3. Let M > 1. The symbol rM{x X) of Rm (cf. (3.3)^ satisfies in local 
coordinates: 

(3.14) \D^D^rM{x,^,X)\ < (sin^)-2^+i-l«l-l^l (O'^-^-l"! /x)"'^. 



Proof. We have that 



The terms in the parametrix symbol ^i>oq-d-l are constructed as solutions to the 
successive equations for m eNq: 



'|Q:|+A;+t=m " ^ I 



1 for m = 0, 

for m = 1, 2, . . . , 



cf. e.g. Seeley [S67], (1). We can use the truncated composition formula in Lemma 3.2 to 
compute the symbol rM of Rm with expansions in up to M homogeneous terms: 

By (3.15), 

Thus Fm consists of the following terms: 

Using the estimates (3.1) and (3.2) together with \D^D^pa-k{x.,^)\ < we 
obtain from Lemma 3.2 with d' = 2d, di = —d — I and d2 = d — k that for Z > 1: 

I D^DfrM-k-l {q-d-h Pd-k)\ 

(3 J6) < (sin ^yM+k+l-2l-\a\-\l3\ ^^^2d-d-l+d-k-{M-k-l)-\a\ ^y2d 

< (sinvp)-2^+i-l«l-l^l(e)''^-^-|'^l(^,/x)-2^ 

since k-l> -M + 1. 

For Z = we find in view of Lemma 3.2 2°, since k < M, 

\D^D^rM-kiq-d,Pd-k)\ 

(3.17) < (sin<^)-l-(^-'=)-|«|-|/3|(^)2d-d+d-fe-(M-fe)-|a|^^^^^-2d 

< (sin ^)-^-i-l«l-l^l(0'''-'^-|"l(^, . 
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Finally, 



(3.18) 



< (sin^)-2^+i-l«l-l^l(0^-'^-l"l(e,;u)-'^ . 



An addition of the contributions (using that ^ 1) gives (3.14). □ 

Now let us show how the kernel of V^{t) is estimated. Similarly to Corollary 2.6, it 
will be convenient to write P = — sUq, where Ho is the orthogonal projection onto 
the zero eigenspace of P, and £ > is chosen < the lowest positive eigenvalue, whereby 
P^ = P + eUo is > e. Here Hq is the ipdo of order with kernel Yl'^^i^j{x)(fj{y)*., for an 
orthonormal basis (fi, . . . of the zero eigenspace. Then 

V{t) = V^t) + (1 - e-^*)no, where V^{t) = e"*-^'; 

and it is the latter operator that needs investigation. V^{t) is defined from the resolvent 
Ql^Qx-{e- A)-ino = {P' - A)-i by 

(3.19) V%t) = ^ [ e-'^QldX. 



JC 

When t G C+ with argument argt = 6 E [0, say, we must assure that Re At — )■ — oo on 
the integral curve C. With ^ —0 denoted 2(fi, this is assured if A runs on a contour formed 
of the rays A = re^*''', connected near by a circle of radius e' < e passing to the right of 
0. Denote infx^c Re A = ci > 0. 

For simplicity of notation we drop the £-index in the next calculations. 

Remark 3.4. The argument 9 = f — of t runs in [0, f [ , when (p runs in ]0, ^j. On this 
interval, sm(p = sin(2</?) = cos^, so they can be used interchangeably in our estimates. 

To find Vlj{t) we can plug Q'j^ into the integral (3.19) and study the kernels, trying 
to get an estimate in terms of \t\. This cannot quite be achieved with (3.19), but better 
estimates are obtained if we first make use of the resolvent formula 



Qx = -X-' + X-^QxP. 



This gives 
(3.20) 



V{t) = ^ f e-'\-X-'+X-^QxP)dX=^ [ e-'^X-^QxPdX, 
Je Jc 

dtV{t)^-^ J^e-'^QxPdX. 



Thus dtJCy^ is the kernel of the integral of the M-th remainder —{Q\P)']^ of —QxP- We 
know that ICv' vanishes at t = 0, and want to show boundedness of the last integral 
applied to the kernel of the M-th remainder. Here (cf. also (3.3)) 

QxP={y2Q-d-i + QM)P = y2, Q-d-iP + RMQxP = y2, O^d-iP + RMPQx- 

KM 
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The term RmPQx already belongs to the M-th remainder {Q\P)'m^ ^^id its kernel is 
estimated by 

II^M^'||£(L2,/f"+i)(sin¥'|A|)-^ 

where RmP can be treated by another application of Lemma 3.2. 

We shall here need an estimate of L2-bounds in terms of symbol seminorms. Many 
variants are known, and we use the following, found in Marschall [M87], Theorem 2.1. 

Theorem 3.5. Let a e 5°o(^" x that for some C > 0, TV e Nq with iV > f , 

and all a, ^ e N[J with <'|a| < < |^| < I, 

(3.21) sup,,^(0-l"l|i^fi?fa(x,OI < C < oo. 

Then the associated operator A = Op(o) is bounded on L2{W^), and || A||£(L2(M")) ^ ^■ 
The dependence of the operator norm on C follows from an inspection of the proof. 

Proposition 3.6. The symbol r = ro(rM,p) of R = RmP satisfies 

(3.22) \D^D^r\ < (sin^)-2^+i-l«l-l'3|(^)2d-A^-|a|(^^^)-d. 
Moreover, for M > n + 1 + 2d, 

^3 \\R\\ciL,,Hr^+^) < (sin(^)-2^-3"/2-2(A)-i, and hence 

\\RMPQx\\c(L,,Hr^^^) < (sin<^)-^^-3"/2-3(A)-i|A|-^ 



Proof. An application of Lemma 3.2 1° to the composition R — RmP shows (3.22). Next, 

= (E|„<„+Ji^^° «"«.,))* = (E|„<„+Ji«'iii(.,))*. 

where an application of Lemma 3.2 3° shows that for |7| < n + 1, 

\D^D^r^x,^,X)\ < (sin(^)-2A^+i-("+i)-l/3|-l«l(0'''-^+("+')-l"l(A)-^ 

According to Theorem 3.5, R'^ is bounded in L2 when M > 2ci + n + 1, with operator 
norm estimated by the seminorms of the symbol derivatives with \a\ < N, \^\ < 1 for some 
N > n/2. Hence the norm is estimated by 

(sinyp)-2^-"-i-("/2+'5)(A)-i < (sin(^)-2^-3V2-2^A)-i, 

where d equals ^ or 1. Summing over I7I < n + 1, we obtain the first line in (3.23), and 
the second follows, since ||(5a||£(L2) — (sin(/? lAI)"-*^. □ 

For the term Ylii<MQ-d-iP we can apply Lemma 3.2 to the compositions, much as in 
Proposition 3.3. 
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Proposition 3.7. For the M-th remainder R'^ = {^i^MQ-d-iP)'M ofYl,i<MQ-d-iP, 
the symbol of^^^j^Q-d-lP is in local coordinates a sum of homogeneous terms of degree 
0, —1, . . . , —M + 1 plus a remainder r'j^ of the form 

(3.24) t'm = q-dP'u + ^m, 
where 

(3.25) rM(a;,e,A)| < (sin^)-2^+i-l"l-l^l(02'^-^-|"l(e,/i)"^^ 

Proof. We shall use estimates from the proof of Proposition Y. Recall from (3.13) that p 
only differs from p in the principal term. We have 

q-d-l opk-d = , , , , , D^q^d-ld^Pd-k + rM-k-l{q-d-hPd-k), 

^—^k+l+\a\<M ^ 

where the terms in the sum are homogeneous of degrees > — M, and the remainders are 
estimated as in (3.16) for I > 0, (3.17) for I = 0, contributing to (3.25). It remains to 
consider 

(3-26) J2i^M^-d-l °p'm = I];<^^o(9-d-/,PM)- 

For I > 0, ro{q-d-hp'M) satisfies estimates as in (3.25). For I = 0, 

ro{q-d,p'M) = q-dPM + riiq-d,p'M): 

where ri{q-diP'M) satisfies estimates as in (3.25) by Lemma 3.2 2°. We collect the terms 
in 

then (3.24) holds with (3.25). □ 

Then we can finally show the estimate of the remainder kernel: 

Theorem 3.8. Let P he selfadjoint strongly elliptic of order d > on M , with 7(-P) > 0. 
The remainder kernel K,v^ satisfies for M > n+l+2d, t e C+ withoxgt = f— 2<^ e [0, 

(3.26) |/Cv^(^,y,t)| < (sin<^)-^^-3-/2-3e-'^<^^min{|t|,l}, 
where c' > if-f{P) > 0, c' = if-f{P) = 0. 

Proof. If 7(-P) > we use the preceding estimates directly to analyse 

Here {QxP)'m = RmPQx + R'mi Propositions Z and V. They contribute as follows: 
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Prom RmPQx we get using (3.23) 
(3.28) 

< (sin(^)-2^-3-/2-3( /'"e-^^°^l*l'-(r)-V-idr+ / e-'''^''\e)-^\e\-^ dip') 

<(sin^)-2^-3-/2-3g-c'Ret^ 

for some c' > 0. 

From R'jyj we get two terms. One is the operator with symbol tm, which satisfies the 
estimates (3.25). Since the kernel equals J^^^z M\z=x-yi we find for M > n + 1 + 2(i that 

|/Cg (a;,y,A)|<(sinv.)-2^+i(A)-^ 

Rm 

Insertion in the integral as in (3.28) gives a bound (sin(/?)~^-^"'"^e~'^ 
To find the contribution from q_dP'Mj ^^^^ perform the integration: 

^^e-*^(/-A)-VMdA = e-*^VW- 

Here we can borrow an estimate from Lemma 3.9 below. By (3.33), 

\D^D^{e-'p'p'M)\ < (sin(/;)-^-l«l-l'3|(^)-^-|«le-^^^*<?)' 

< (sin^)-^-l«l-l^l(0-^-|«le-'^'^^^ hence 
I^Op(e-.°VM)l ^ (sin¥')"'^e-'='^^* when M > n+ 1. 

Since the latter estimates are dominated by that from RmPQx, we conclude that 
\dtlCv' {x,y,t)\ is estimated as in (3.28). Then an integration with respect to t using 
that ^i/^(x, y, 0) = shows (3.26). 

In the case where 7(-P) = 0, the above considerations will be valid for V^{t) as in (3.19). 
We then have to add (1 — e~'^*)no, which has a smooth kernel bounded by min{|t|, 1} and 
we reach the conclusion in the theorem. □ 

Now we turn to the contribution from the homogeneous terms. 

Lemma 3.9. Let M let ai, . . . , aM be nonnegative integers with 

(3.29) a = ai + --- + aM>l, 

and let f{x, ^, A) be a (matrix-formed) symbol of the form 

(3.30) fix, ^, A) = f,{pd - X)-^'f2{Pd - A)-^^ ■■■{Pd- A)-^-/m+i, 

where the fj{x,^) are il^do symbols of order Sj e homogeneous for |^| > 1. Denote 
si + ■ ■ ■ + sm+1 = s, then the order of f is k = s — ad. Let Fx = Op{f{x,^,X)) on W^, 
and let E{t) be the operator family defined from Fx for Ret > by 

(3.31) E{t) = ^ I e-'^FxdX, 
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when t = e**^2 Then E{t) = Op(e(a;, t, ^)), where the symbol 

(3.32) e{x, t,^) = ^J^ e-'^f{x, C, A) dX 

satisfies: 

(i) e(x,s-'^t,sO = s"+'e(x,t,0/or|e|>l, s>l, 

(3 33) 

(ii) |Dfi:>^e(x,t,OI < (sm^)-'^-|"|-|^l(0^+''"'"'e-"^^'<^>'. 
T/ie kernel of E{t) satisfies for d + k > —n 

(3.34) \)CEix,y,t)\ < (sin(^)-^-('^+'=+")/'^e-'^'^^*|tr('^+'^+^)/^ 

with c' > 0. 
If(T>2, 

(3.35) |DfDfe(a;,t,e)| < (sin(^)-^-l"l-l^l|t|(0^'^+''"l''le-'^'^<'^^^>' . 
/n t/izs case, the kernel satisfies for d + k < —n 

^, • . ^ ,j,,(\t\(\logRet\ + l)ifd + k = -n, 

Proof. As in [G96], Lemma 4.2.3, we can pass the operator definition through the integral. 
To estimate e, we first consider |^| < 1. We use the residue theorem and that is 
self adjoint to obtain 

CT— cRet 



lit, X, e) I = I ^ ^ e-'^f{x, e, A) d\\ < (1 + |r-i)e 



Here, c = -i{p^{x,^)). 

For 1^1 > 1, we replace C by a closed, homogeneous curve Cc^c around the spectrum of 
^). Cc,c coincides with C on a annulus of inner radius c\^\'^ and outer radius Cl^l*^, is 
closed by the segments of the boundary of this annulus which lie to the right of C. Then 
by homogeneity 



Ht,xM = [ e-*V(^,e,A)rfA| < (sin<^)--(e)'^(0'e-i^^*l«l 



Combining the two estimates, we conclude (3.33) for a = ^ = 0. The derivatives 
D^D^e{x, t, ^) are of a similar form, with k and a replaced by A; — \a\ resp. d + |q;| + \/3\. 

To show (3.34) foTd + k> —n, we estimate K-e by comparing e with its homogeneous 
extension e'*: 



lCE{x,y,t)= [ e'^'=-y>^e^{x,t,^)d^+ [ e^(^-^>^(e - e'^) 
Jk" J\£.\<i 
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Using (3.33) and a homogeneous variant, 

\JCE{x,y,t)\ < (sin^)-'^e-'=i^^* / g"'^^ (f^ 

+ (sinv?)-"e-'^i^^* / e-'^^^^*l^l'((0^+^ + 1^1'^+'')^?^ • 

The first integral is = (Ret)-('^+'=+")/'^ = (sin(^)-('^+'=+")/'^|t|-('^+'=+")/'^, while the second 
remains bounded as \t\ — )■ 0. 

Now consider the case where a > 2. As |/| < (A)~^ away from M_|_, the integral converges 
uniformly in t > 0. We may deform C to a closed curve in the left half-plane, where / is 
holomorphic, to conclude e{x,0,^) = 0. Also, using {—X){pd — A)~^ = 1 — Pd{Pd — A)~^, 

dte{x, t,O^^J^ e-'^ {-X)f{x, A) dX 

can be expressed in terms of e and a second term of the same form, with one of the Sj 
replaced by Sj + d. By (3.33) 

\dte{x,t,^)\ < (sin(^)-'^(e)^'^+'=e-^^^*<^>' 
and hence, since the value at i = is 0, 

|e(x,t,OI < (siny?)-"|t|(0^"'+''e-'^^^*<^>' . 

This shows (3.35) ioi a — (5 — Q. The proof for D^D^e{x,t,^) is analogous. The estimate 
(3.36) is obtained similarly to (3.34), using (3.35) instead of (3.33). □ 

Theorem 3.10. 1° In local coordinates the kernel terms satisfy for some d > 0; 

{(sin(/9)(^-'^)/^|t|('-^)/^ ifd-l>-n, 
|t| (|logRet| + 1) zfd-l = -n, 
\t\ if d — l < —n. 
2° Let N — max{^, ^ + Ad + 7} , Then the full kernel satisfies 

\}Cv{x,y,t)\ < sin((/?)-^e-^'^'^*|tr"/^ . 

For a given cq > we can modify to satisfy inf^;^^ 7(^^(0;, ^)) > cq/ then c' can be any 
number in ]0, co[ . 

Proof. 1°. For / > 1, the assertion follows from Lemma 3.9, (3.34) resp. (3.36). 
On the other hand, we explicitly compute for / = 

\K:v_^{x, y,t)\ = {27ry\ f e^^^-^'^^e-*^"^^'^) 

<g-ciRet^ /" g-C2Ret|pO(a:,OI _^ f (^£-02 Ret|p°(a;,6l _ g-C2 Ret|p';i(x,0| ) 

Jk" •^I^I<i 
<g-ciRet(^ f e-R<^*l?l'd$ + l) < e-^'^^*(ReO""/^ • 

The assertion follows, since Ret = |t| simp. 

2°. We choose M = n + 1 + 2d in Theorem 3.8 and add }Cv_a-i ioi <l < M. The 
most singular terms dominate. 

See Section 2 for how to obtain the allowed range of c'. □ 
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Theorem 3.11. 1° In local coordinates, K-V-d satisfies for some c' > 0; 

(3.38) \]Cv_a{x,y,t)\ < (sin^)-^-^e-'^' ^''^Itl |a; - yp'^-". 
For I > 1, the kernels )Cv_d-i satisfy 

(3.39) \}Cv_,_Ax,y,t)\< 

\t\ \x — y^^-d--"^ if d — I > — n, 
(sin^)-["-^+^]+-2'e-'=' ^^*<J \t\ (I log\x-y\\ + l)ifd-l = -n, 



\t\ if d — I < —n. 



2° Let N = ^ + 4:d + 7. Then the full kernel satisfies 

\ICv{x,y,t)\ < sin(^)-^e-'='^^^|t| \x - y\-'^-'' . 

If ^(P) > 0, we modify p'^ to satisfy inf^;,^ 7(p°(a:, ^)) > 7(-P), then d can he any number 
m]0,7(P)[. 

Proof. 1°. For Z > 1, we obtain from Lemma 3.9, (3.35), that 

\D^V-d-i{x,t,0\ < (sin(^)-2'-l"l(0''-'-l"l|t|e-'='^^* . 

The estimate (3.39) then follows from the kernel estimates in Proposition 2.2. 

For / = 0, v_d{x,t,^) = e-^P and we obtain (3.38) as in the proof of Theorem 2.4. 

2°. We choose M = n + 1 + 2d in Theorem 3.8 and add ICv_a-i for < / < M. The 
most singular terms dominate. □ 

The estimates hold also when (/? e [0, ^[ is replaced by (/? e ] — ^, 0] and sint/? is replaced 
by |sinv3|. Note that for t with argt = ^ g] — f,f [, cos^ = |sin2(/p| = |sin(^|, when 
d = ±(^ — 2\(p\) (of. Remark 3.4). Then we can formulate the final result obtained by 
combining Theorem 3.10 and Theorem 3.11 as in the proof of Theorem 2.5, as follows: 

Theorem 3.12. Let P be selfadjoint strongly elliptic of order d > on M , with ^(P) > 0. 
The heat kernel JCy satisfies for t e C+ with axgt = $ & ] — ^, ^[ the Poisson estimate, 
where N = max{^, ^ +4d + 7}: 

(3.40) \ICvix,y.t)\ < (cos^)-^e-^(^)^-^ (d(a;,y)+ |t|i/rf)<^ + '^''^'^"^ + 

Proof. In the region where < d(x, y), d{x^ y) = d{x^ y) + Itl"^/*^, and the asserted esti- 

mate follows from Theorem 3.11. In the region where > d{x, y), we use \t\^/'^ = d{x,y) + 
and Theorem 3.10. □ 

Remark 3.13. The +7 in N partly stems from repeated rounding up to the nearest 
integer. It may be reduced by choosing M + 2d closer to n and using a version of Theorem 
3.5 with a non-integer number of derivatives, or sharper versions for negative orders. 

Our result applies in particular to the Dirichlet-to- Neumann operator. For this operator 
ter Elst and Ouhabaz [E013] have estimates in terms of — A'^-th powers of cos 6, where the 
dimension n enters nonlinear ly in N . 
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4. Kernels of heat semigroups for perturbations of fractional 
Laplacians and the Dirichlet-to-Neumann operator 

This section complements the general upper bounds from Section 2 with lower estimates 
in the case of fractional powers of the Laplacian and the Dirichlet-to-Neumann operator. 

Let A be the (nonnegative) Laplace-Beltrami operator on the closed, compact Riemann- 
ian n-dimensional manifold M ; it defines a selfadjoint nonnegative operator on L2(M), also 
denoted A. In this case, is an elliptic pseudodifferential operator of order d on M, 

with positive principal symbol |^|^, defining a selfadjoint nonnegative operator on L2(M); 
it generates a holomorphic semigroup = e~*^'^^^ with C°°-kernel for t > 0, 

lCya{x,y,t) = {5,,V\t)8y). 

The semigroups e"*"^ and V'^{t) are related by subordination formulas. For d — 1^ they 
assume a simple form: 

Lemma 4.1. Let A > 0. One has for t > 0; 

1 2 

Proof. Let a = tVX/2 and let x = |s~^; then dx = —^s~ids, and equation (4.1) is 
turned into 

/"CO 2 2 

(4.2) = / e"^ -S^2cZx. 



To show this, note that the left-hand side I{a) satisfies I{a) e C^(M_|_), limQ,^>o+ (ck) = 
^/^T, and for a > (with y = ax~^, dy — —ax~'^dx): 

dal{a) = / e-"" -^{-Aa)x-^ dx = -2 I e~^~^ 2dy=-2I{a). 
Jo Jo 

Thus I ( a) = ce-^'" with c = ^/7^ . □ 

By Zolotarev [Z86] (see also Grigor'yan [G03]), there exists for any < d < 2 a non- 
negative function r]f{s) such that 

nOO 

(4.3) e-'^"^' = / e-^^ Vti^)ds- 

Jo 

Here r]f has the following properties 

(4.4) r/f(.) = t-2/V(^) {s,t>0), 

(4.5) r]f{s) <ts-^-i {s,t>0), 

(4.6) r]^{s)^ts-^-i (s>t2/d>o). 
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By an application of the spectral theorem, we obtain for alH > 0, 

/•OO 

(4.7) V\t)f = e-*^''V = / e--^/ vf{r) dr , for all / e H^{M). 

Jo 

In view of (4.7), it holds that 

poo pCXl 

{5,,V\t)5y) = {5,, / e-^'^Sy rjfir) dr) = / (4, e'^^) rjfir) dr, 
Jo Jo 

resulting in an identity for the kernels: For alH > 0, 

poo 

(4.8) ]Cv^ {x, y,t)= / /Ce-rA {x, y) 4{t) dr , for {x, y) e M x M. 

Jo 

Using this formula, we can deduce upper and lower estimates for ICyd from those known for 
/Cg-TA. The following upper and lower estimates are well-known (see e.g. L. Saloff-Coste 
[SlOl); 



(4.9) ^^e-^.^<K.-,.(x.y)<^^^e-^«^' 

V[X,y/T) V[X,y/T) 

Here V(a;,r) denotes the volume of a ball of radius r around x. For a closed compact 
n-dimensional manifold M, V{x,r) = r"' for small r, and V{x,r) equals the volume of the 
connected component containing x when r > diamM. Hence 

(4.10) V(x,V^)-i = (W2)-i + l. 

Theorem 4.2. Let < d <2. The kernel of the semigroup V^it) = e-*^'^' satisfies for 
t > 0; 

(4.11) /C^_,../. {X, y) = ^^^^J^^r,ay v) + ^'Z')- + l) • 

Proof. The upper estimate follows already from Corollary 1.6. The following proof more- 
over extends to give the lower estimate. Inserting the heat kernel bounds (4.9), (4.10) into 
(4.8) and using (4.5), we find 

ICy4x,y,t)< / {T-^/^ + l)r^^{T)e-^'^dT 
Jo 

<t T T ^ e ^ -r dr + t T^2e*^ - dr . 
Jo Jo 

By a change of variables r i->- Cd{x, y^r, the first term equals 
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Similarly, the second term is 

r-OO 



/•OO 

t{Cd{x,yf)-i T-i-^ e-^l^ dr = 
Jo 

and altogether, 



t 



d{x,y) 



d 



t 

d{x,yy 

On the other hand, using the uniform bound /Cg— ta (x, y) < r""/^ + 1 and (4.4), we obtain 



nOO 

= / (t-"/rfr-"/2 + 1) r/f (r) dr = t"^/'^ + 1 . 
Jo 



Thus 

t 



/Cv^.(a;,y,t)< min{t-/'^ + l,-— — + . 



d{x,yy 

If t^/'^ > d{x,y), 



^-n/rf<^-n/d(^^^ + l^ " = t{d{x,y) + t'/'')-^-'^ 



and 



On the other hand, for t^^^ < d{x, y) we have d{x^ y) = d{x^ y) + 1^/'^ and hence 

^-^{dix, y)-" + 1) < — / ( idix, y) + tVd)-n + _ 

d{x,yy^ ^ >- {d{x,y) + t^/d)dV ^ ' ) 

This shows " < " in (4.11). 

To show the opposite inequality in (4.11), note that the integrand in (4.8) is non- 
negative, and (4.9), (4.10) imply 

lCya{x,y,t)= lC,-r.{x,y)rjf{T)dT> / {r-^/^ + 1) rjfir) e-""^ dr 

Jo Ja 



for a = max{t2/'^,rf(a;,y)2}. Now, for r > d{x,yf, ^-c^^^ > ^-c rj^j^^^^ ^^_g^^ 

/•OO 

}Cvd{x,y,t)> / (r-^'/^ + l) tr^-^ dT = t(a~'^ +a"^) 
= minjt-'^/^, td(a;, 2/)"""'^} + min{l, td{x, yy^} 

> t{d{x, y) + t^/dyn-d ^ ^^^^^^ ^ ^1/dyd q 
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For d = 1 this complies well with the explicit kernel formula (2.3) for the Poisson 
operator solving the Dirichlet problem for the Laplacian on Rlj:"''"'^. 

We also consider the case where M is the boundary of a compact (n+l)-dimensional Ric- 
mannian manifold M with boundary. With A denoting the nonnegative Laplace-Beltrami 
operator on M, we shall compare /C _tv^ with the kernel of the semigroup generated by 
the (nonnegative) Dirichlet-to- Neumann operator Pdn on M. Pdn is the operator map- 
ping u to the normal derivative dyU, where u is the harmonic function on M with boundary 
value u. It is known (cf. [G71]) that Pdn is an elliptic pseudodifferential operator of order 
1 on M with the same principal symbol as . 

bmce is a classical strongly elliptic ijjdo of order d, Theorem 2.5 applies to all 

operators of the form P = A*^/^ + P' with P' classical of order (i — 1, giving upper 
estimates of the absolute value of the kernels; note that no selfadjointness is required. For 
such operators we can also show lower estimates. 

Theorem 4.3. Let d e ]0, 2[ and let P be a classical i(^do of order d with the same principal 

symbol as A^/^. Then the kernel ofV{t) = e~*^ satisfies for t > 0; 

(4.13) 

\K:v{x, y,t)\<t {{dix, y) + t^^) — ^ + y) + t^/dyd^ + ^-'^^H {d{x, y) + t^/df-n-d^ 

for any c\ < ^{P) (ci = ^{P) if Corollary 2.6 applies). Moreover, there is an r > such 
that 

(4.14) \ICv{x, y,t)\>t {d{x, y) + t^/dyd-r^^ for d{x, y) + 1^''' < r. 

Proof. As P and A'^/'^ have the same principal symbol, 

V{t) = V^it) + V, 

where V is of lower order, more precisely V is the difference between the first remainders 
for V{t) = e~*^ and V'^{t) = e~*'^'^^^ , as in the second line of (2.25). Hence 

(4.15) \]Cv'{x,y,t)\<e-'''H{d{x,y) + t^/''f-''-''. 

Now (4.11) and (4.15) together imply (4.13). 

To obtain the lower estimate (4.14), we note that 

(4.16) cs""-'^ - c's^-"-'^ = cs-"-''(l - c'c'h) > 2-^cs-"-^ when s < c/(2c'), 

so for t in a bounded set where e~'^^* < c', the lower estimate in (4.11) implies that (4.14) 
holds for small d{x, y) + t^/*^. □ 

We can also obtain upper and lower estimates for the Dirichlet-to-Neumann operator. 

Theorem 4.4. The kernel o/e~*^^^ satisfies for t > 0; 

(4.17) /Ce-ti-^^ (x, y, t) < ^^^^ {{d{x, y) + t)"" + l) , 

and there is an r > such that it satisfies 

(4.18) JC^-tPuN (x, y,t)>t {d{x, y) + t)"'"", for d{x, y)+t<r. 

Proof. Here Pdn is known to be selfadjoint nonnegative, and the semigroup has real, 
nonnegative kernel ([AM07], [AM12]), so that we may omit absolute values. The upper 
estimate (4.17) follows from Corollary 2.6. The lower estimate (4.18) follows from Theorem 
4.3 since Pdn differs from A^/^ by a classical ipdo of order 0. □ 
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Remark 4.5. This work was inspired from a conversation of the second author with W. 
Arendt and A. ter Elst in August 2012, where we suggested the apphcabihty of pseudo- 
differential methods as in [G96] to the Dirichlet-to-Neumann semigroup. We have very 
recently learned of the efforts of ter Elst and Ouhabaz in [E013], giving an analysis of the 
Dirichlet-to-Neumann semigroup by somewhat different methods, and obtaining some of 
the same results as those presented here. 
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